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Preface 


The  objective  of  this  thesis  was  to  determine  plasma 
electron  number  density  and  electron-neutral  collision 
frequency  from  the  measurement  of  the  complex  Impedance  of 
a  radio  frequency  (RF)  plasma  discharge  within  a  barrel 
type  reactor  under  specified  operating  conditions  (gas  type, 
pressure,  RF  power  level). 

This  report  is  limited  to  a  specific  plasma  model  and 
impedance  measurement  technique.  The  choice  of  model  and 
experimental  configuration  was  found  to  limit  the  accuracy 
of  Impedance  measurements  v^en  applied  to  or der-of -magnitude 
ranges  of  frequency  or  power.  Application  to  the  general 
class  of  inductively  driven  plasma  reactors  shows  the  most 
promise  of  success  when  the  specific  reactor  can  be  isolated 
from  stray  capacitive  effects. 

I  would  like  to  express  appreciation  to  Dr.  Peter 
Bletzinger  and  the  Air  Force  Wright  Aeronautical  Labs  for 
suggesting  and  sponsoring  this  study  and  to  Lt  Col  William  F. 
Bailey,  for  his  professional  insight  and  assistance.  I  would 
also  like  to  thank  Us.  Barbara  A.  Graham  and  my  brother, 
William  F.  Hermanson,  for  their  graphical,  clerical,  and 
moral  support. 


(This  document  was  typed  by  Us.  Barbara  A.  Graham) 
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Abstract 


I 

The  plasma  within  an  RF  driven  inductive  coil  inter¬ 
acts  with  the  electrical  generation  circuit  providing  feed¬ 
back  useful  for  control  of  bulk  plasma  parameters.  An 
investigation,  employing  a  tunable  RF  circuit,  was  made 
applying  complex  impedance  measurements  as  a  plasma  diag¬ 
nostic.  Plasma-circuit  effects  are  modeled  as  capacitive 
coupling  to  th<‘  driving  inductor;  the  magnitude'  of  the 
plasma  model  capacitance  being  proportional  to  a  compli'x 
plasma  permittivity.  Associated  Impedance,  when  annulled 
by  tuning,  results  in  a  real  resistive  1  f)ad  depende-nt  upon 
plasma  characteristics.  Impedance  diagnostics  avoid  Lhr' 
use  of  traditional  internal  probes.  Application  of  this 
technique  to  inductive  reactors  is  found  to  be  inff'rior 


to  capacitive  reactors. 


I.  Introduction 


Background 

Plasma  reactors  are  used  extensively  in  the  manufacture 
of  integrated  circuits  and  thin  film  solar  collectors.  With¬ 
in  th(?  reactor  a  plasma  is  formed  when  energy  is  transferred 
to  a  low  pressure  gas  through  the  application  of  a  radio 
frequency  (RF)  electric  field.  Subsequent  molecular  dis¬ 
sociation,  ionization  and  particle  redistribution  make 
possible  the  etching  or  deposition  of  desired  substances 
onto  receiver  substrates  used  for  integrated  circuit  pro¬ 
duction  . 

Plasma  reactors  applicable  to  RF  processes  are  distin¬ 
guished  by  physical  geometry  and  method  of  excitation.  Two 
reactor  configurations  which  are  used  in  industry  are  the 
parallel  plat*'  type  and  the  barrel  type.  The  latter  lends 
itself  to  application  in  large  scale  production  processes. 
Figure  1  is  a  schematic  of  two  reactor  types  used  by  the 
Ih'wlett-Packard  Corporation.  Reactant  gas  is  bled  into  a 
chamber  being  evacuated  by  pumps.  An  RF  varying  potential 
is  applied  to  conducting  plates  causing  an  active  plasma 
region  to  form.  The  positions  of  substrates  or  wafers  are 
shown  on  which  the  desired  etching  or  deposition  process 
can  be  aff<*ct<'d.  The  geometry  of  the  barrel  type  offers 
an  obvious  advantage  of  being  volume  efficient.  Here  the 
wafers  ca  be  stacked  such  as  dominoes  and  the  reactor 


Figure  1.  Plasma  Reactors 


designed  to  accommodate  a  large  number  of  wafers.  A  high 
degree  of  uniformity  results  from  the  use  of  such  reactors 
(Ilic,  D.B.,  1982:20).  The  method  of  RF  excitation  may  be 
accomplished  through  a  capacitive  or  inductive  coupling. 

The  latter  finds  particular  use  with  the  barrel  type  reactor. 
Figure  2  shows  schemati*cs  of  these  alternate  excitation 
methods.  The  Inductively  coupled  reactor  displays  a  large 
plasma  region  unobstructed  by  internal  electrodes. 

Characterization  of  reactor  parameters  in  relation  to 
the  electrical  driving  values  of  frequency  and  power  provides 
feedback  useful  in  the  control  of  plasma  processes.  Compli¬ 
cated  molecular  gases  such  as  Silane  (SiH^)  are  often  used 
because  of  their  varied  chemical  reactions,  of  interest  to 
integrated  circuit  technology.  The  level  and  frequency  of 
the  applied  power  strongly  influences  the  energy  of  the 
electrons  and  this  selectively  determines  the  dissociation 
products  produced  as  well  as  the  deposition  or  etch  rate. 

The  method  of  RF  excitation,  similar  to  driving  an 
antenna,  requires  that  the  impedance  of  the  load  be  matched 
to  that  of  the  circuit  in  order  to  transfer  the  maximum 
amount  of  power  to  the  load.  RF  plasma  reactors  used  in 
both  research  and  industry  have,  therefore,  impedance 
matching  networks  incorporated  into  the  driving  circuits  to 
perform  this  function.  It  is  possible  to  capitalize  on  this 
feature  to  make  plasma  measurements.  Since  the  values  of 
plasma  parameters  influence  the  impedance  matching  circuit 
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design.  It  is  impossible  to  give  an  all-purpose  matching 
network  for  a  general  plasma  (Norstrom,  H.,  1979:349).  In 
a  recent  paper,  H.  Norstrom  explains  four  possible  classes 
of  matching  networks  and  selection  criteria.  Although  the 
actual  circuits  may  vary,  the  optimization  of  power  transfer 
through  impedance  matching  allows  identification  of  a  criti¬ 
cal  point  at  which  measurements  may  be  taken. 

Impedance  is  an  electrical  circuit  parameter.  Similar 
to  resistance,  its  use  arises  from  the  fact  that  in  an  AC 
circuit  the  voltage  and  current  vary  periodically  and  one 
may  lead  or  lag  the  other  in  phase.  The  circuit  impedance 
is  a  complex  quantity  consisting  of  both  magnitude  and  phase 
Information  and  may  be  expressed  in  polar  form  or  converted 
to  rectangular  form  to  express  separately  a  real  and  ima¬ 
ginary  part.  The  circuit  impedance  can  be  directly  measured 
using  equipment  appropriated  for  this  experiment.  Alterna¬ 
tively,  individual  circuit  components  can  be  calibrated, 
thereby  avoiding  repetitive  measurements.  Impedance  conveys 
information  concerning  the  circuit  load,  in  this  case,  the 
plasma.  By  Judicious  measurement,  impedance  may  be  used  as 
a  plasma  diagnostic.  This  has  been  demonstrated  in  a  recent 
article  (Ilic,  D.B.  1981)  where  impedance  techniques  were 
developed  and  applied  to  a  capacitive,  parallel-plate  RF 
plasma.  In  the  study  which  follows,  similar  techniques  are 
developed  and  applied  to  an  inductively  driven  barrel-type 
reactor.  Since  the  impedance  measurements  are  performed 
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on  existing  circuitry,  they  have  the  advantage  of  being 
external  to  the  plasma  avoiding  the  use  of  traditional 
probes  or  microwave  techniques  and  hence,  are  nondisturbing 
to  particle  motions  or  internal  electric  field  distributions. 

Problem  Statement 

M  ■  I  —  II  ■  ■  « 

The  problem  investigated  in  this  study  demonstrates 
that  the  presence  of  an  inductive  RF  plasma  within  a  barrel- 
type  reactor  interacts  with  the  electrical  driving  circuit. 
This  interaction  is  identified  in  the  form  of  the  complex 
Impedance  of  a  circuit  element  and  can  be  used  as  a  diagnos¬ 
tic  to  determine  plasma  electron  number  density  and  electron- 
neutral  collision  frequency  for  specified  operating  conditions 
(gas  type,  pressure,  RF  power  level). 

Scope  and  Assumptions 

Study  was  initially  conducted  using  a  Helium  plasma, 
chosen  because  of  its  relatively  simple  atomic  structure, 
well  documented  collisional  cross  sections  and  it  character¬ 
istically  uniform  behavior  over  wide  ranges  of  pressure  and 
energy.  This  provided  a  fundamental  reference  for  develop¬ 
ment  of  models  and  evaluation  of  the  diagnostic  techniques. 
Time  constraints  did  not  permit  a  desired  expansion  of  this 
study  to  include  molecular  gases  such  as  Silane  (SIH^)  which 
are  commonly  used  in  industrial  Silicon  disposition  processes. 

The  pressure  range  considered  was  bounded  by  0.01  torr 
at  the  lower  limit  due  to  the  pumping  capacity  of  the  vacuum 
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system  roughing  pump.  An  upper  bound  of  1.0  torr  was  select¬ 
ed  because  of  reported  poor  quality  of  thin  films  deposited 
at  such  pressures  and  attendant  high  power  requirements. 
(lliCy  D.B.,  1982:5).  The  plasma-generating  electric  field 
is  determined  by  the  power  and  frequency  of  the  driving  RF 
signal.  The  signal  generator  was  operated  over  the  range 
from  10  to  35  MHz  and  the  amplifier  permitted  continuous 
operation  at  power  levels  up  to  40  watts  average  power. 

The  spatial  and  temporal  nature  of  the  electric  field 
within  the  plasma  was  not  studied  in  detail  but  approaches 
that  of  a  cylindrical  winding  with  negligible  end  effects. 
This  can  be  visualized  as  a  uniform  axial  magnetic  field 
and  an  azmuthal  electric  field,  both  varying  sinusoidally  in 
time.  These  approximations  do  not  distort  results,  since 
Impedance  values  are  calculated  from  the  dielectric  response 
of  the  plasma.  The  dielectric  response  is  a  measure  of  the 
Intoto  response  of  plasma  constituent  particles  under  the 
influence  of  the  applied  field  resulting  in  a  net  induced 
molecular  dipole  moment  (Marion,  J.B.  ,  1980:23).  Readers 
Interested  in  a  more  rigorous  electromagnetic  field  analysis 
are  referred  to  Tang  and  bis  study,  applicable  to  RF  fields, 
using  emissive  probes  (Tang,  J.T.,  1979:1458). 

Two  necessary  models  are  incorporated  in  this  study  to 
relate  the  plasma  to  an  electrical  circuit  measurement  (the 
Impedance).  The  first  is  formed  by  treating  the  plasma  as 
a  lossy  dielectric  which  assigns  to  the  plasma  a  complex 
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permittivity  dependent  upon  plasma  parameters.  This  is  a 
well  accepted  treatment,  following  from  electromagnetic 
field  theory  using  Maxwell's  equations  to  develop  an  expli¬ 
cit  expression  for  a  complex  propagation  constant.  A  par¬ 
ticularly  thorough  analysis  is  given  in  a  text  written  by 
Marr  (Marr,  G.V.,  1968:21).  The  second  model  arises  from 
circuit  analysis  wh<»n  the  plasma  is  considered  as  an 
isotropic  m(‘dium  contained  within  a  circuit  element. 

Genera  1  Approach  and  Organization 

In  the  section  to  follow,  general  theory  is  presented 
to  explain  plasma  phenomenon  and  establish  a  background  for 
developing  th(?  mathematical  relations  necessary  to  define 
and  quantify  electrical  circuit  values  and  plasma  parameters 
Specific  plasma  parameters,  electron  number  density,  and 
electron-neutral  collision  frequency,  are  selected  as  most 
descriptive  variables.  These  are  investigated  to  determine 
their  functional  dependence  upon  the  RF  driving  frequency 
and  host  gas  pn*ssure. 

Th«>  sp<*ciri<-  equipment  used  and  th('  physical  configu¬ 
ration  art!  described  and  shown  in  Figure  3.  A  glass  barrel- 
typt!  plasma  tub«!  was  used  with  an  associated  low  pressure 
vacuum  system  and  gas  supply.  An  impedance  matching  elec¬ 
trical  network  was  incorporated,  using  an  RF  tunable  signal 
generator  and  linear  amplifier,  coupled  to  the  plasma  tube 
with  an  inductive  coil.  A  monochrometer  and  digital  photome 


ter  was  available  permitting  spectral  analysis  of  the  emitted 
plasma  radiation. 

Procedures  are  developed  that  yield  values  for  the 
plasma  parameters  from  corresponding  circuit  impedance 
measurements.  During  operation »  a  frequency  is  selected,  a 
plasma  is  generated,  and  the  electrical  matching  network 
tuned  to  match  the  load  impedance  to  that  of  the  driving 
generator.  Pressure  and  power  levels  are  then  varied  while 
the  standing  wave  ratio  (SWR)  meter  is  monitored  to  deter¬ 
mine  the  Impedance  matching  condition , 'that  is,  when 
maximum  power  is  being  transferred  to  the  plasma  and  a 
minimum  is  being  reflected  back  through  the  electrical 
network . 


>1. 


Figure  3.  Equipment  Layout 
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I I .  Theory 


This  study  considers  an  RF  generated,  weakly  ionized, 
non-equilibrium  plasma.  Prior  to  the  application  of  the 
electric  field,  the  gas  contained  within  the  reactor  exists 
in  a  nearly  equilibrium  state.  The  effect  of  an  imposed 
field  is  to  cause  existing  free  electrons  (and  ions),  re¬ 
sponding  to  a  Lorentz  force,  to  gain  kinetic  energy.  The 
presence  of  neutral  particles,  unaccelerated  by  thf'  field, 
inhibits  electron  motion  and  facilitates  energy  transfer 
among  particles.  As  the  electrons  attain  energies  suffi¬ 
cient  to  overcome  ionization  thresholds  during  collisions, 
a  plasma  is  generated.  For  any  given  level  of  driving 
power,  a  steady  state  is  reached  where  the  energy  and 
particle  losses  from  diffusion,  recombination  and  radiation 
balance  the  energy  being  absorbed  by  the  plasma  from  the 
field.  Implicitly,  particle  collisions  perform  a  control¬ 
ling  function  in  both  the  electron  production  and  loss 
mechanisms  occurring  within  the  plasma. 

The  effect  of  an  alternating  electric  field  on  a  free 
electron  is  to  cause  the  electron  to  oscillate  within  limits 
of  velocity  and  kinetic  energy  determined  by  the  field  strength 
and  period.  Over  time  the  electron  will  not  gain  energy  from 
the  field  under  these  conditions.  In  the  presence  of  other 
particles,  however,  resulting  collisions  redirect  the  motions, 
allowing  the  electrons  to  accumulate  random,  thermal  energy. 


If  the  applied  electric  field  is  sinusoidally  varying  in 
time  and  expressed  as:  E(t)  =  E  cos  ut ,  the  rate  (energy 
is  gained  from  the  field  can  be  approximated  by: 


•  ^  =  V 

dt  ~  2m  (u)''  +  ' 

(Stamm,  M.R. ,  1981:26) 

where  e  and  m  are  the  electron  (charge  and  mass  resp«'c t  i  vel y , 
E  is  the  field  amplitude,  w,  the  driving  rre(|uency  and  v , 
the  collision  frequency. 

The  collision  frequency  is  important  because  it  also 
appears  in  energy  and  particle  loss  terms.  The  rate  of 
free  electron  production  is  directly  related  to  the  rate  of 
energy  shown  and  thus  dependent  explicitly  on  the  collision 
frequency.  Two  limiting  cases  may  be  considered.  In  the 
DC  approximation  w  approaches  zero  and  the  rate  of  produc¬ 
tion  scales  inversely  with  the  collision  frequency.  In 
this  instance  the  mechanism  may  be  thought  of  as  an  electron 
being  accelerated  by  the  field,  suffering  a  collision  re¬ 
sulting  in  the  exchange  of  its  energy,  then  being  reaccelc'r- 
ated  until  subsequent  collisions.  The  distance  hc'tween 
collisions  is  termed  the  mean  free  path  and  the  rat('  of 
collisions,  the  collision  frequency.  Clearly,  the  occur¬ 
rence  of  collisions  acts  to  impede  the  progress  of  th<' 
electrons  to  some  limiting  velocity,  termed  the  drift 


velocity. 

In  the  case  where  the  electric  field  is  rapidly  oscil¬ 
lating,  w  is  much  greater  than  the  collision  frequency  and 
the  rale  of  production  depends  directly  upon  the  rate  of 
collisions  but  is  strongly  inversely  dependant  upon  the 
driving  frequency.  Each  of  these  limits  must  be  considered 
in  this  (jxperiment  since  the  operating  ranges  of  control 
parameters  (gas  pressure,  power,  frequency)  place  the 
driving  frequency-collision  frequency  relationship  within 
the  transition  between  these  two  regions.  In  the  absence 
of  compensating  effects,  experimental  data  can  be  expected 
to  reflect  the  effects  of  this  transition  between  adjoining 
regimes . 

Because  of  their  more  rapid  response  to  the  applied 
field,  it  is  the  electrons  rather  than  the  ions  which  drive 
plasma  proc(.'sses.  The  rate  at  which  energetic  electrons 
stimulate  particular  reactions  is  best  contemplated  by 
assigning  to  c!uch  collisional  event  a  specific  cross  section 
which  implicitly  relates  the  probability  of  occurrence  to 
a  1  unction  of  the  average  energy  of  the  electron.  For 
mo  I  t'c.u  I  ar  gases,  energy  dependant  electron  impact  cross 
sections  are  oft»>n  v(?ry  complicated  functions.  Fortunately, 
the  chcjice  of  lh?lium  for  a  constituent  gas,  provides  a  rela¬ 
tively  complete  set  of  collision  cross  sections.  Electron 
excitation  and  ionis'.ation  cross  sections  for  Helium  have 
been  accurately  measured  as  functions  of  electron  energy 
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and  comprehensive  data  is  readily  available  therefore  it  is 
not  necessary  to  repeat  here.  (Brown,  S.C. ,  1966:133) 

Electrons  below  inelastic  energy  thresholds  suffer  only 
elastic  collisions  with  the  molecules.  As  the  electrons 
gain  kinetic  energy  and  excitation/ionization  thresholds 
are  reached,  inelastic  collisions  occur  causing  radiative 
losses  through  excitation/de-excitation  events  and  the 
electron  population  increases  through  ionizations.  As  the 
el(ictron  number  density  grows,  a  characteristic  plasma  glow 
is  observed.  Since  the  available  power  level  is  limited, 
balancing  loss  mechanisms  limit  fractional  ionization  levels 
in  this  experiment  to  the  order  of  10”^, 

From  this  discussion,  it  is  apparent  that  the  rate  of 
plasma  generation  is  determined  by  the  electron  number 
density  and  the  frequency  of  collisions.  The  collision 
fr«?quency  is  expressed  as  the  electron  velocity  divided  by 
the  im'an  frt-e  path.  This  can  be  related  to  the  cross  section 
as : 


V  =  vno  (2) 

(Chen  ,  F.F.  ,  1977:137) 

where  v  is  the  t'leetron  velocity,  n  is  the  neutral  particle 
number  density  and  o  is  the  electron  impact  cross  section. 
The  neutral  particle  number  density  is  a  function  of  host 
gas  pressure.  The  mean  electron  velocity  can  be  expressed 
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in  terms  of  energy  and  the  collision  frequency  then  written 


V  -  BpoTgy*  (3) 

(Hirsh,  N.H.,  1928:177) 

% 

where  T  is  energy  in  electron  volts,  p  is  gas  pressure 

©V 

in  torr,  o  is  the  cross  section  and  B  is  a  proportionality 
constant  equal  to  21.0  x  10^^  for  these  units. 

For  select  atomic  gases,  including  Helium,  the  colli¬ 
sion  frequency  is  found  to  be  relatively  constant  over  a 
broad  range  of  energy.  Hydrogen  also  belongs  in  this  group 
but  its  use  was  considered  too  volatile.  Heg,  a  mixture 
of  Helium  and  mercury  belonging  to  this  category  was  un¬ 
available.  Brown  has  compiled  measurements  of  these  gases 
and  gives  the  result  for  Helium  as: 

V  -  2.37  X  10’  p  (4) 

(Brown,  S.C.,  1966:168) 

Electron  losses  are  also  controlled  by  the  collision 
frequency.  Random  motion  caused  by  the  thermal izing  effect 
of  collisions  bring  particles  to  the  container  walls  where 
they  recombine.  In  addition  to  this  "free  diffusion",  the 
collective  behavior  of  the  plasma  introduces  ambipolar 
diffusion  in  that  the  higher  mass  of  the  ions  and  attendant 
slower  response  Inhibit  the  electron  diffusion.  If  the 
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electrons  and  ions  were  to  diffuse,  each  at  its  own  free 
diffusion  rate,  the  net  charge  difference  would  continue 
to  grow  over  time.  Instead,  the  plasma  responds  by  estab¬ 
lishing  internal  fields  so  as  to  reduce  the  electron  flux 
and  increase  the  ion  flux  until  both  species  diffuse  at  the 
ambipolur  rate.  Ambipolar  diffusion  dominates  in  plasmas 
that  have  physical  dimensions  that  are  large  compared  to  the 
Debye  length.  In  this  experiment,  the  Debye  length  is  on 
the  order  of  a  millimeter,  much  less  than  the  tube  diameter 
of  10  cm. 

The  effect  of  ambipolar  diffusion  is  to  increase  the 
ion  diffusion  rate  by  a  factor  approximating  the  ratio  of 
the  ion  to  electron  energy  (Stamm,  M.R.,  1982:35).  We  will 
consider  a  fac;tor  of  40  to  be  appropriate  here.  The  co¬ 
efficient  of  diffusion  is  approximated  by  the  following 
express  ion : 


AT  ’ 

I.  vA  „  3  ev 

D  =  ^-  =  b  X  10’ 


(Stamm,  M.H.,  1982:48) 


(5) 


where  A  is  the  atomic  weight  in  atomic  mass  units  and  the 
remaining  (iuantiti(?s  are  as  previously  defined.  For  a  1  torr 

O  *  1  if  ^  ^ 

Helium  plasma  at  300  K  with  a  cross  section  of  5  x  10  cm  , 
the  ambipolar  diffusion  coefficient  is  calculated  to  be 
100  cm^/sec  (S'lmm,  M.R. ,  1982:50). 
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This  is  in  agreement  with  experimental  data  (Brown,  S.C 
1967:128).  The  time  for  diffusion  can  be  written  in  terms 
of  a  characteristic  diffusion  length: 


1  .  D 
T  “  XT 

(Stamm,  H.R.,  1982:44) 


(6) 


The  diffusion  length  is  based  on  the  geometry  of  plasma 
container  and  is  found  by  solving  Laplace's  equation  for 
the  steady  state  case.  For  a  cylindrical  geometry  it  is 
found  to  be: 


1 

TJ 


+  ( 


2.405 

R 


2 

) 


(Hirsh,  M.N.,  1978:188) 


(7) 


Applying  the  above  to  this  experiment  gives  a  value  for  A 
of  2.07  and  a  time  for  diffusion  of  0.04  sec. 

The  relatively  slow  response  of  the  ions  introduces 
another  important  effect  which  impacts  the  ability  of  the 
electric  field  to  sustain  the  plasma  and  causes  frequency 
dependant  behavior  in  the  plasma  parameters.  Initial  dis¬ 
placement  from  an  electron-ion  equilibrium  position  results 
in  a  restoring  force.  Particle  inertia  causes  restoring 
particle  motions  to  overshoot  the  equilibrium  positions, 
giving  rise  to  an  oscillation  which  occurs  at  a  character- 
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istic  frequency.  By  expressing  the  particle  motions  as  an 
equation  of  motion  and  applying  Maxwell's  equations,  the 
plasma  frequency,  is  found  to  be: 


»_ 


4yne- 

m 


(Chen,  F.F.,  1977:73) 


(8) 


Noting  the  dependence  on  the  electron  number  density,  it  is 
convenient  to  combine  the  constants  and  use  the  following 
approximation : 


fp  -  9000  n^  (9) 

(Chen,  F.F.,  1977:73) 


Maxwell's  equations  also  leads  to  a  dispersion  relation 
for  transverse  electromagnetic  waves  propagating  through  a 
plasma.  Considering  a  simple  case  where  there  are  no  ex¬ 
ternal  magnetic  fields  imposed,  the  dispersion  relation  can 
be  written  in  terms  of  a  propagation  constant,  K,  and  the 
velocity  of  light,  C,  as: 

K*  -  ^  (o)*  -  u  *)  (10) 

C‘  p 

(Chen,  F.F. ,  1977:102) 


The  salient  feature  of  this  relation  is  that  the 
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propagation  constant  becomes  imaginary  when  the  wave  fre¬ 
quency  exceeds  the  plasma  frequency.  Stated  another  way; 
when  the  plasma  frequency  is  fixed,  due  to  a  specified 
number  density,  electromagnetic  waves  below  the  plasma 
frequency  will  not  propagate  through  the  plasma.  For  a 
density  of  10^”  cm"’,  the  plasma  frequency  is  on  the  order 
of  a  gigahertz.  This  explains  the  wide  use  of  microwave 
radiation  for  plasma  study. 

At  frequencies  below  the  plasma  frequency  the  extent 
of  penetration  of  the  electromagnetic  radiation  is  limited. 

A  measure  of  the  depth  of  penetration  is  termed  the  skin 

• . 

depth,  6,  which  describes  the  distance  required  for  the 
radiation  to  be  attenuated  to  1/e  of  its  Initial  value  and  " 
obeys  the  following  relation: 


(Chen,  F.F.,  1977:103) 

Here  C  is  the  velocity  of  light  and  other  terms  are  as 
previously  defined.  This  indicates  that  at  all  frequencies 
below  the  plasma  frequency  there  is  finite  penetration. 

This  discussion  presents  a  paradox  when  applied  to  the 
RF  case  of  study.  It  would  appear  that  the  RF  energy  being 
used  to  generate  the  plasma  would  not  be  allowed  to  pene¬ 
trate  the  boundry  once  the  plasma  was  formed  to  sustain  the 


discharge.  Although  this  holds  true  for  the  proper  boundry 
conditions,  in  this  experiment,  the  dimensions  of  the  plasma 
(10  cm)  are  small  compared  to  the  free  space  wavelength  (15 
meters),  (Sanders,  S.G.,  1978:2690).  With  this  consideration 
and  allowing  for  a  finite  skin  depth,  sufficient  penetration 
is  made  to  sustain  the  discharge  (Uarion,  J.B. ,  1980:151). 

This  discussion  again  points  out  the  importance  of 
the  particle  density.  This  parameter  acts  to  balance  the 
free  electron-ion  production  caused  by  the  field  with  the 
losses  due  to  diffusion  to  and  recombination  at  the  container 
walls  by  controlling  collisions. 

During  operation,  the  plasma  is  formed  and  allowed  to 
stabilize  before  measurements  are  taken.  In  this  steady 
state,  the  rate  of  free  electron  losses  must  equal  the  rate 
of  production.  The  equation  of  continuity  is: 

^  -  IT  “  -  <1*> 

In  the  first  term,  n  is  the  neutral  gas  number  density  and 

ft 

is  the  average  rate  of  ionization  found  by  summing  indiv¬ 
idual  ionization  rates  from  the  ground  and  excited  state 
colllsional  events.  The  second  term  considers  losses  from 
diffusion  as  previously  discussed.  The  third  term  is  inclu¬ 
ded  to  account  for  volume  recombination  events.  Since  this 
term  is  non-linear  with  respect  to  electron  number  density. 


♦ 
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tbe  equation  of  continuity  can  be  solved  in  the  steady  state, 
for  tbe  number  of  free  electrons  in  the  plasma.  This  becomes: 


“a'l  - 


iF 


(13) 


% 

An  Increase  in  driving  power  level  causing  an  increase  in  the 
production  rate  of  ionizations  is  shown  by  the  above  express¬ 
ion  to  imply  an  increase  in  number  density.  This  trend  is 
an  expected  observation  in  tbe  experimental  data. 

It  is  necessary  to  relate  the  plasma  parameters  to  the 
electrical  circuit  in  order  that  electrical  measurements 
will  yield  meaningful  information  about  tbe  collision  fre¬ 
quency  and  number  density.  This  linkage  is  affected  by 
modeling  the  plasma  as  a  lossy  dielectric  with  a  permittivity 
explicitly  described  by  the  plasma  parameters.  The  complex 
permittivity  is  given  by  either  of  the  following  two  ex¬ 
pressions  : 


(Illc,  D.B.,  1981:1544) 


1 


qm 
01^  + 


1  gym 

<i)(a)i  +  V* ) 


(15) 


where  a  is  a  numerical  constant  arising  from  the  dependence 
of  the  plasma  frequency  on  number  density  and  equals  3182 
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In  liKS  units.  Other  variables  are  as  previously  defined. 

The  second  equation  is  in  rationalized  form  delineating  the 
real  and  imaginary  parts  of  the  function. 

The  plasma  permittivity  can  be  incorporated  into  ex¬ 
pressions  for  a  measurable  circuit  parameter,  the  impedance. 

« 

Tbe  impedance  is  also  represented  by  a  complex  nuniber  al¬ 
though  for  purely  capacitive  or  inductive  circuit  components, 
the  expressions  are  pure  imaginary  numbers.  If  the  lossy 
dielectric  plasma  is  considered  to  be  tbe  medium  within  a 
circuit  capacitor,  the  resulting  capacitance  is  found  from 
the  basic  air-gap  capacitance  by  using  the  permittivity  as 
a  proportionality  expression.  The  product  of  the  Imaginary 
portion  of  the  plasma  permittivity  with  the  imaginary  capa¬ 
citive  impedance  produces  an  additional  real  term  analogous 
to  resistive  components. 

The  Impedance  may  be  given  in  polar  form  thereby  carry¬ 
ing  magnitude  and  phase  information.  This  is  a  more  usable 
form  when  impedance  matching  is  considered  since  the  desire 
here  is  to  maintain  continuity  of  phase  between  the  circuit 
and  load.  Polar  form  was  also  compatible  with  the  available 
measurement  equipment. 

Circuit  Model 

The  model  chosen  to  represent  plasma-circuit  relation¬ 
ships  is  dependant  upon  the  reactor  geometry  and  type.  In 
the  simplest  case  of  the  capacitive  type  parallel-plate 
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reactor  investigated  by  111c,  the  plasma  fills  the  volume 
between  the  plates  and  the  capacitance  is  given  by  the 
expression:  C  -  where  is  the  basic  capacitance 

of  the  reactor  determined  geometrically  from  the  area  and 
separation  of  the  parallel  plates.  (Ilic,  D.  B.  1981:1544) 
It  follows  immediately  that  the  reactor  impedance  is  given 
by: 


1 

iwe 

P  o 


(16) 


In  the  case  of  the  inductively  excited  barrel  reactor 
used  in  this  experiment,  the  circuit  model  is  more  involved. 
Here  the  value  of  the  basic  capacitance  cannot  be  determined 
geometrically.  With  this  in  mind,  the  model  chosen  to  repre¬ 
sent  the  plasma-circuit  relationship  is  that  of  a  capacitor 
in  parallel  with  the  driving  inductor.  For  generality,  the 
total  capacitance  associated  with  the  inductor  is  separated 
into  components  and  accounts  for  the  effect  of  parasitic 
capacitance  as  well  as  plasma  effects.  Figure  4  shows 
schematically  the  plasma  model  and  its  relationship  to  the 
driving  circuit.  As  shown,  L  represents  the  inductor,  C^, 
the  portion  of  the  parasitic  capacitance  unaffected  by  the 
plasma  and  C^,  the  capacitance  interior  to  the  plasma  tube 
which  is  modified  by  the  presence  of  a  plasma  of  relative 
permittivity,  e^.  For  this  analysis,  the  plasma-circuit 
model  is  treated  as  an  unknown  impedance,  Z^.  When  condtilned 
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with  the  driving  inductor,  the  circuit  reduces  to  an 
equivalent  impedance,  Z  ,  at  the  point  in  the  circuit 
where  measurements  are  to  be  extracted. 


Figure  4*  Plasma  Uodel 

The  expression  for  the  impedance  of  this  model  is: 


_ 1 

i«(Cg  +  CpC^) 


(17) 


The  analysis  appropriate  to  each  of  the  above  models 
differs  only  slightly,  however  the  parallels  drawn  are 
significant.  The  analysis  of  the  capacitive  case  can  be 
gleaned  from  a  critical  review  of  Ilic's  study.  In  this 
article,  Ilic  presents  a  method  of  determining  electron- 
neutral  collision  frequency  and  electron  number  density  in 
a  capacitive  type  plasma  reactor.  Mr.  Ilic  did  not  address 
an  equally  Important  class  of  plasma  reactors,  that  of  the 
inductive  type.  The  question  of  whether  this  method,  with 
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some  modification,  could  be  extended  for  use  with  an  induc¬ 
tive  type  plasma  reactor  suggested  the  investigation  presented 
in  the  body  of  this  thesis.  Since  motivation  for  much  of  the 
analysis  stems  from  parallel  arguments,  more  than  a  cursory 
sumnary  of  this  reference  is  warranted. 

It  must  be  emphasized  that  the  method  presented  was 
applied  to  capacitive  type  reactors  only  and  indeed  is 
applicable  only  to  a  particular  subset  requiring  a  maximum 
energy  transfer  into  the  plasma.  This  requirement  is  satis¬ 
fied  when  a  SWR  of  1.0  is  established.  The  plasma  reactor 
Is  treated  as  an  electrical  circuit  element  for  analysis. 

The  other  essential  elements  are  the  RF  signal  generator 
and  matching  network.  When  the  matching  network  is  properly 
tuned,  the  imaginary  part  of  the  complex  impedance  vanishes 
and  the  signal  generator  drives  a  known,  purely  resistive 
load. 

i  The  matching  network  consists  of  two  variable  capacitors 

i  - 

and  an  inductor.  The  plasma  load  is  also  capacitive.  Using 

( 

‘  traditional  circuit  analysis,  the  electrical  components  are 

reduced  to  an  expression  for  a  single  complex  impedance  based 
on  the  values  of  the  individual  circuit  elements  which  are 
then  measured  or  determined  from  calibrated  curves. 

Choosing  the  plasma  lossy  dielectric  model  with  permit¬ 
tivity  given  by  Equation  (14)  and  using  Equation  (16),  the 
impedance  can  be  rewritten  as: 
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(18) 


(w*  -  ivw)(vw  -  i(w*  -  Wp*)] 
"  u>6^[(vw)^  +  (u^- 


The  numerator  can  be  expanded  and  simplified  into  the  form: 


Numerator  *  +  iw(u  *  -  u*  -  v*)]  (19) 

P  P 


Similarly,  the  denominator  gives  the  form: 


<i)  2 

Denominator  =  - D*  + 


(20) 


Writing  the  impedance  in  rectangular  form  gives  a  real  and 
imaginary  component  such  that  Zp  =  Rp  +  iXp.  Also,  a  factor 
G  is  defined  in  the  following  manner: 


O  -  <  ^  -  D'  ♦  i 


(O' 


(21) 


This  mathematically  arbitrary  but  convenient  form  allows 
simplified  expressions  for  the  impedance.  The  components  are 
as  follows: 


R_ 


V(0p* 

(o'Ci^G 
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(22) 
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These  components  yield  an  Impedance  from  values  of 
plasma  parameters.  Our  desire  is  to  do  the  inverse »  that  is» 
knowing  the  value  of  Impedance,  discern  the  value  of  plasma 
quantities.  We  note  that  the  equations  are  coupled  and 
therefore  solving  for  explicit  expressions  of  plasma  para¬ 
meters  is  tedious.  Again  an  arbitrary  choice  is  made  and 
a  new  quantity,  H,  is  defined  by: 


Careful  examination  of  the  expression  R^,  solved  for  v, 
reveals  that  this  factor  H  must  also  obey  the  following: 


Although  not  immediately  apparent,  several  intuitive  substi¬ 
tutions  can  show  that  the  above  two  expressions  lead  to  an 
identity.  In  the  process  of  performing  this  critical  step, 
an  explicit  expression  for  collision  frequency  arises. 

Further  manipulation  yields  an  expression  for  plasma  fre¬ 
quency  which  we  recall  is  dependant  upon  number  density.  The 
final  expressions  for  these  parameters,  in  agreement  with 
Ilic's  study,  are  the  following: 
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(26) 
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The  experimental  aesthetics  of  this  method  arise  from 
the  realization  that  once  calibration  curves  for  the  circuit 

t 

components  are  determined,  plasma  collision  frequency  and 
number  density  can  be  calculated  from  the  matching  network 
counter  readings.  This  is  particularly  convenient  for  the 
configuration  presented  in  this  article  which  embodied  an 
automatic  tuning  network. 

Mr.  Ilic  presented  one  example  to  verify  that  reasonable 
values  could  be  obtained  for  a  Nitrogen  discharge  driven  at 
a  frequency  of  13.56  MHz.  Lack  of  a  more  dynamic  range  en¬ 
courages  further  investigation.  Variable  components  were 
measured  at  low  frequency  which  makes  wider  application 
somewhat  suspect  unless  considerable  attention  was  given  to 
eliminating  stray  capacitive  effects.  Even  small  unaccounted- 
for  values  of  inductive-capacitive  coupling  introduce  signi¬ 
ficant  Increases  in  apparent  component  values  as  studied  by 
R.  A.  Olsen  and  observed  by  this  thesis  student  (Olsen,  R.A., 
1961:2).  Collision  frequency  is  determined  from  a  resistive 
term  which  arises  from  the  presence  of  the  imaginary  term 
in  the  expression  for  the  plasma  permittivity.  The  magnitude 
of  this  resistive  term  is  small  and  demands  high  measurement 
accuracy  to  infer  an  associated  collision  frequency. 

Returning  to  the  case  under  Investigation  in  this  study. 
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a  parallel  development  can  be  performed.  Using  the  same  ex¬ 
pression  for  permittivity  and  the  expression  for  the  impedance 
applicable  to  the  inductive  barrel -type  reactor,  the  combi¬ 
nation  gives  the  following: 


(u*  -  ivfa)) 


(28) 


For  convenience,  a  numerical  constant  y  is  defined  as: 


Y 


(29) 


By  redefining  the  term  G  used  previously,  the  impedance  can 
again  be  expressed  in  component  form.  Let: 


I 

(30) 


It  is  instructive  to  compare  the  above  expression  to  that 
pertaining  to  the  capacitive  case.  Equation  (21)  and  note 
the  inclusion  of  the  constant,  y.  Likewise,  the  term  appears 
in  the  following  expressions: 


VU)  * 


(31) 
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(32) 


Again  we  wish  to  solve  for  the  plasma  parameters  and  with 
the  same  Justification  as  before,  we  define  a  new  H,  such 
that : 


H  «  yR 


P 


irXp’  +  -e- 


oiC. 


(33) 


This  selection  gives,  after  some  careful  manipulation,  the 
sought  after  expressions  which  explicitly  relate  the  colli¬ 
sion  frequency  and  number  density  to  measured  impedance 
values.  These  relations  are: 


(34) 


(35) 


In  the  above,  a  is  a  constant  arising  from  the  conversion 
from  plasma  frequency  to  collision  frequency  {a  =  3.182  x 
10*  cm*/sec*).  It  remains  to  deter  ine  numerical  values 
for  the  real  and  imaginary  parts  of  the  plasma  impedance. 
This  is  accomplished  through  measurements  of  the  inductor, 
L,  and  the  equivalent  impedance,  Z  .  Circuit  analysis  re- 
duces  Zp  to  the  following  expression: 
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(36) 


-  (dLX  +  iojLR 

V  =  _ _ ^ 

p  -  R  +  i((i)L-X  ) 
e  e 

A  vulu<’  Tor  L  was  found  to  be  2.70  microhenrys.  The  ex- 
prt'Ksion  for  Z  can  be  rationalized  to  give  the  real  and 

I-  p  r. 

iinat^inary  parts,  however,  the  present  form  is  usable  in  a 
suhi-oul  inc  of  a  computer  program  which  is  used  for  data 
reduction.  (St'e  Appendix  C). 
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III.  Equipment 


The  equipment  employed  in  the  study  is  discussed  under 
the  following  functional  categories :  Plasma  tube,  vacuum 
system,  electrical  system,  optical  system,  and  an  overall 
configuration.  • 

Plasma  Tube 

The  plasma  tube  diagrammed  in  Figure  5  is  a  locally 
fabricated,  low  pressure,  pyrex  glass  tube  having  a  diameter 
of  10  cm  and  a  length  of  56  cm.  The  tube  incorporates  oppo¬ 
sing  input-output  openings  to  provide  a  non-turbulent  gas 
flow,  a  pressure  port,  and  a  viewing  tube  and  window.  Two 
similar  plasma  tubes  were  used.  One  without  a  viewing  tube 
was  used  for  study  of  atomic  gases;  the  second,  includes  the 
viewing  tube  provision  which  is  necessary  to  prevent  contam¬ 
ination  of  the  viewing  window  by  RF  deposition  when  molecu¬ 
lar  gases  such  as  Silane  are  studied. 
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Vacuum  System 

The  vacuum  system  utilizes  a  DuoSeal,  Model  1397,  low- 
vacuum  mechanical  pump.  Positive  shutoff  valves,  controllable 
flow  valves,  mixing  manifold  and  associated  flow  meters  and 
pressure  gauges  control  the  flow  and  content  of  the  gas  con- 
tained  within  the  plasma  tube.  Provisions  are  made  to  dilute 
the  vacuum  pump  exhaust  thereby  avoiding  possible  volatile 
reactions  which  can  occur  between  ambient  air  and  molecular 
gases  such  as  Silane.  Since  this  experiment  uses  a  flowing 
gas  system,  the  roughing  capacity  of  the  vacuum  pump  estab¬ 
lishes  a  practical  lower  pressure  limit  of  0.01  torr.  The 
vacuum  system  is  shown  schematically  in  Figure  6. 

The  gas  under  investigation  flows  from  bulk  gas  supply 
cylinders  through  a  pressure  regulator  which  is  adjusted  to 
maintain  a  pressure  consonant  with  downstream  meters  and 
valves  (approximately  10  psi).  The  gas  passes  through  a 
calibrated  chamber  where  a  transducer  is  used  to  determine 
an  accurate  volume  flow  rate.  The  gas  then  passes  through 
a  glass  flow  tube  providing  a  secondary  flow  rate  reading. 
Controllable  flow  valves  and  positive  shutoff  valves  admit 
the  desired  flow  into  a  manifold  which  stabilizes  the  flow 
into  the  plasma  tube.  A  direct  reading  pressure  gage,  read¬ 
able  to  one  torr,  provides  approximate  system  pressure.  A 
pressure  sensor  calibrated  for  the  specific  gas  under  study 
is  connected  close  to  the  active  plasma  region  to  give 
accurate  readings  down  to  10""  torr.  Exhaust  is  expelled 
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through  a  constant  volume  vacuum  pump  and  diluted  to  ambient. 
Electrical  System 

The  electrical  system  makes  use  of  a  Hewlett  Packard, 
Model  606A,  RF  signal  generator,  a  Hewlett  Packard,  Model 
5245L,  digital  frequency  counter,  a  Heathkit,  Model  2040 
antenna  tuner,  a  linear  amplifier,  and  associated  meters  to 
deposit  RF  energy  into  the  plasma  by  means  of  an  induction 
coil.  A  component  demanding  stringent  specifications  for 
this  experiment  is  the  linear  amplifier.  The  Amplifier 
Research,  Model  40L,  linear  amplifier  employed  here  incor¬ 
porates  unique  circuitry  which  allows  operation  into  any 
load  Impedance.  This  is  critical,  since  plasma  generation 
Introduces  an  unknown  impedance  and  a  prolonged  mismatched 
condition  often  exists  until  the  tuning  procedure  can  be 
accomplished  and  measurements  taken. 

An  HP  4815  RF  Impedance,  Vector  Voltmeter  was  used  for 
calibration  and  impedance  measurements.  Figure  7  shows  the 
electrical  system  components.  The  tuner  is  shown  expanded 
and  labeled:  "matching  network",  because  of  its  importance 
in  circuit  analysis. 

Optical  Diagnostic  System 

A  SPEX  Model  1870  Czerny-Turner  Spectrograph  provided 
a  means  of  correlating  emission  spectra  to  calculated  plasma 
parameters.  The  spectrometer  with  associated  variable  scan¬ 
ning  drive,  digital  photometer,  photomultiplier  and  cooling 
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Flffuire  7«  Electrical  System 


system,  and  strip  chart  recorder  is  coupled  to  the  plasma 
tube  with  a  choice  of  optical  doublets.  The  spectrometer 
has  an  f/l  of  8.6,  entrance  and  exit  slit  widths  variable 
from  0  to  3  mm  and  a  choice  of  three  slit  heights.  Because 
of  limited  table  space,  an  appropriate  input  doublet  was 
procured  corresponding  to  each  slit  height  to  focus  the 
plasma  emission  onto  the  entrance  slit  without  the  need  for 
changing  overall  object  to  image  distance.  In  each  case, 
the  input  optics  would  slightly  overfill  the  slit  to  provide 
maximum  throughput  with  minimum  background  interference. 
Figure  8  shows  the  input  optics  relationships. 

Configuration 

The  overall  equipment  configuration  is  shown  in 
Figure  9.  The  plasma  tube  complete  with  mixing  manifold 
was  mounted  on  a  platform  designed  for  ease  of  alignment 
relative  to  the  spectrometer.  The  electrical  system  is 
coupled  to  the  plasma  tube  via  the  induction  coil.  The 
optical  system  is  coupled  through  the  input  optics.  The 
plasma  tube  and  driving  coil  were  enclosed  within  a  grounded 
aluminum  cage.  External  coaxial  electrical  connections 
were  added  to  facilitate  Impedance  measurements. 
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Plasaa  Tube 


Figure  8.  (^tical  System 


Figure  9.  Equipment  Layout 


IV.  Procedure 


Impedance  measurements  were  taken  using  the  following 
procedure.  With  the  vacuum  pump  maintaining  a  constant- 
volume  pumping  rate,  the  flow  control  valve  was  adjusted  to 
establish  a  desired  pressure  within  the  plasma  tube.  The 
frequency  of  interest  was  selected  on  the  signal  generator 
and  RF  impedance  meter.  Both  were  calibrated  using  the 
digital  frequency  meter  to  assure  equivalence.  The  power 
level  of  the  amplifier  was  raised  to  produce  a  plasma.  If 
a  plasma  was  not  readily  generated  due  to  a  highly  mis¬ 
matched  condition,  the  matching  network  was  readjusted. 

Once  a  plasma  was  generated,  and  an  approximate  forward 
power  level  established,  the  Standing  Wave  Ratio  (SWR) 
meter  was  monitored  while  the  matching  network  was  tuned  to 
a  SWR  of  1.0. 

Ibe  actual  matching  procedure  is  dependant  upon  indiv¬ 
idual  load  parameters.  The  tuner  has  three  variable  compo¬ 
nents  but  no  consistent  order  of  adjustment  was  universally 
productive.  However,  familiarity  reduced  tuning  time  con¬ 
siderably.  In  tuning  an  unknown  impedance,  the  most  reliable 
procedure  was  to  reduce  the  sensitivity  of  the  SWR  meter  to 
keep  an  on-scale  reading,  adjust  the  tuner  inductor  until 
a  slight  dip  was  observed  in  the  SWR,  then  position  both 
capacitors  for  a  minimum  SWR  and  then,  alternately  increment 
the  Inductor  and  the  antenna  matching  capacitor  while  checking 
for  a  lower  minimum  SWR  reading  after  each  pair  of  adjustments. 
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Once  a  SWR  of  1.0  was  attained,  a  measurement  could  be  taken 


With  the  amplifier  positioned  to  standby,  the  coaxial 
cable  was  removed  from  its  connector  at  the  driving  coil. 

The  impedance  of  the  aggregate  of  forward  components  was 
measured  at  this  point  and  labeled  ^ .  Since  the  Im- 

pedance  was  matched  (SWR  =  1.0),  the  equivalent  impedance 
of  the  coll  and  plasma,  Z  ,  adds  vectorally  with  the  value 
of  to  produce  a  zero  phase  condition.  Thus,  if  we 

measure  i®  polar  form,  Z^  must  necessarily  be  of  the 

same  magnitude  and  opposite  phase.  This  value  of  Z  is  then 

G 

used  in  appropriate  theoretical  equations  to  determine  col¬ 
lision  frequency  and  number  density.  The  following  figure 
indicates  this  measurement  schematically. 

The  value  of  the  driving  inductor,  L.,  was  determined 
from  a  resonance  condition.  The  inductor  has  an  associated 
Interturn  capacitance,  C,  and  the  combination  can  be  treated 


as  a  capacitor-inductor  parallel  circuit  having  an  impedance 


Figure  10.  Impedance  Measurement 
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given  in  magnitude  by: 


|Z| 


(i)L 

l-w^LC 


(37) 


By  making  an  Impedance  ipeasurement  at  several  different  fre¬ 
quencies,  the  value  of  L  was  found  to  be  2.70  microhenrys  and 
C  was  found  to  be  9.5  picofarods. 

As  discussed  in  Section  II,  the  capacitance  measured 
above  was  considered  to  consist  of  an  internal  and  external 
component.  The  sensitivity  of  the  measurements  requires 
that  the  capacitance  be  accurately  apportioned  between  its 
components.  This  was  accomplished  by  taking  impedance 
measurements  with  each  of  two  different  dielectric  materials 
filling  the  plasma  tube.  The  first  material  was  dry  air 
having  a  dielectric  constant  of  unity.  The  second  material 
was  carbon  tetrachloride  which  was  chosen  for  its  dielectric 
constant  of  2.27.  This  value  is  essentially  invarient  to 
frequency  changes  over  the  RF  range.  Since  the  presence  of 
the  dielectric  influences  the  internal  capacitance  only, 
these  measurements  allow  for  solving  simultaneously  for  the 
values  of  the  individual  capacitances.  The  internal  capa¬ 
citance,  C^,  was  found  to  be  0.52  pf  and  the  external 
capacitance,  Cg,  8.98  pf. 
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V.  Results 


Plasma -Circuit  Interaction 

Tlie  thrust  of  this  thesis  is  to  infer  the  magnitudes 
of  plasma  parameters  from  their  effect  upon  electrical 
circuit  elements,  in  particular,  the  driving  inductor.  The 
fact  that  the  presence  of  a  plasma  does  have  detectable 
impact  was  demonstrated  conclusively.  By  generating  a 
plasma  and  tuning  the  matching  network,  a  SWR  of  1.0  was 
attained,  indicating  a  matched  impedance  load.  After  extin¬ 
guishing  the  plasma  (by  momentarily  interrupting  the  power 
at  below  threshold  levels),  a  substantial  deflection  was 
noted  on  the  SWR  meter.  This  indication  of  a  mismatched 
condition,  the  plasma's  effect,  encouraged  attempts  to 
quantify  the  plasma-circuit  perturbation  through  impedance 
measurements. 

Model  Validation 

The  basic  model  chosen  for  analysis  is  that  of  a 
theoretical  plasma-filled  capacitor  connected,  in  parallel, 
to  the  driving  inductor.  This  model  was  shown  to  have  some 
validity  through  the  following  experiment. 

A  variable  capacitor  was  installed  in  a  parallel  posi¬ 
tion  to  the  driving  coil  to  allow  simulation  of  the  plasma. 
When  the  plasma  was  extinguished,  corresponding  to  a  change 
in  the  theoretical  plasma  filled  capacitor,  the  SWR  meter 
deflected  substantially  from  its  annulled  reading. 


The  impedance  matched  conditon  could  be  restored  by  making 
a  small  adjustment  (approximately  2  picofareds)  on  the 
variable  capacitor.  Thus,  changes  in  a  real,  physical 
capacitor,  in  parallel,  duplicate  observed  phenomena  attri¬ 
buted  to  plasma  parameters.  The  converse  suggests  that  the 
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same  effects  can  be  attributed  to  a  theoretical,  plasma- 
filled  capacitor  in  the  same  parallel  configuration. 

Application 

Measurements  were  made  of  the  circuit  impedance  as  a 
function  of  pressure  and  power  using  procedures  outlined  in 
Section  IV.  The  data  was  reduced  using  equations  presented 
in  Section  II  and  Appendices  B  and  C.  The  resulting  plots 
were  then  examined  for  possible  trends  which  would  correlate 
to  the  theory  discussed  in  Section  II.  Figures  11  through 
14  present  examples  of  the  measurements  taken.  The  quantity 
being  measured  is  complex  and  the  meter  readings  are  given 
in  polar  form.  Each  measurement  gives  two  points  of  infor¬ 
mation,  a  phase  angle  and  a  magnitude.  Plotting  the  data 
in  polar  form  was  found  to  be  awkward  for  interpretation 
therefore,  each  measurement  was  converted  to  rectangular 
form.  This  allows  the  real  and  imaginary  parts  of  each 
measurement  to  be  plotted  individually.  Although  the  data 
is  plotted  on  separate  graphs,  the  real  and  imaginary  portions 
must  be  understood  to  be  coupled  in  order  to  maintain  proper 
perspective. 
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Figure  12.  Imaginary  Impedance  Component 


Figure  13.  Real  Impedance  Component 


46 


Examples  of  the  data  show  that  the  impedance  measurably 
responds  to  changes  in  system  parameters  (power «  frequency, 
pressure,  etc.)*  Figures  11  and  12  show  the  impedance 
measured  at  20  MHz.  The  real  portion  exhibits  a  common 
value  at  low  pressure,  and  has  little  dependence  on  pressure 
until  reaching  a  breakpoint  after  which  it  scales  inversely 
with  pressure  to  the  1/3  power.  The  imaginary  portion  shows 
no  dependence  on  pressure  within  measurement  accuracy  limits. 
Figures  13  and  14  show  impedance  measured  at  constant  power. 
Changes  in  frequency  produce  changes  in  both  the  magnitude 
and  slope  of  the  impedance  vs  pressure  curves. 

Figures  15  and  16  show  plots  of  the  admittance  (the 
reciprocal  of  impedance)  as  a  function  of  frequency.  Plotted 
for  comparison  is  the  case  of  an  ideal  resistor-inductor- 
capacitor  (RLC)  circuit  using  constant  component  values 
comparable  to  those  measured.  The  presence  of  the  plasma 
has  a  pronounced  effect,  increasing  with  frequency.  Figure 
16  shows  the  imaginary  portion.  Here  the  ideal  circuit  would 
reach  a  resonance  condition  near  30  MHz.  The  plasma  results 
in  a  resonance  near  22  MHz.  Treating  the  plasma  as  a  lossy 
dielectric  should  produce  the  opposite  effect.  Since  the 
resonant  frequency  is  given  by  u*  =  (LC)“^,  a  reduction  in 
the  capacitance  would  result  in  a  higher  resonant  frequency. 
The  resistance  term  (Figure  15)  arises  in  the  model  from  the 
product  of  the  imaginary  term  of  the  permittivity  with  the 
imaginary  capacitive  reactance.  Measurements  show  a  deviation 


47 


Figure  IS.  Real  Admittance  Component 


from  the  ideal  circuit.  Together  these  results  suggest  the 
plasma-circuit  model  to  be  incomplete. 

Figures  17  through  20  show  reduced  data.  Under  certain 
sets  of  conditions,  the  results  agree  closely  with  theory. 

For  example,  the  slope  of  the  collision  frequency  vs.  pressure 
curve  (Figure  17)  is  expected  by  Equation  (4)  to  be  2.4  x  10*. 
At  20  MHz,  this  value  is  correctly  calculated  at  0.1  torr. 

The  slope  of  the  curve  is  slightly  below  expected  (1.2  x  10*) 
at  lower  pressures  and  slightly  higher  (4.1  x  10*)  at  higher 
pressures  but  the  average  slope  agrees  with  theory  over  the 
entire  pressure  range  investigated.  Similar  results  appear 
for  22  MHz  and  25  MHz.  Other  selected  frequencies,  for 
example,  18  MHz  and  28  MHz  exhibit  a  much  smaller  slope, 
significantly  deviating  from  theory.  Calculated  electron 
number  densities  exhibited  similar  trends.  However,  expected 
uniform  increases  with  respect  to  increasing  power  levels 
were  not  observed.  Nevertheless,  an  increase  in  plasma 
intensity  is  visably  apparent.  To  quantify  this  observation, 
several  Helium  emission  lines  were  spectroscopically  scanned. 

A  plot  of  photoconductor  current  vs.  RF  input  power  is  shown 
in  Figure  21.  Each  line  examined  exhibited  a  continuous 
Increase  over  the  range  of  input  power. 
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Figure  18.  Number  Density 


A0N3n03aJ  NOISmOD 


Figure  20.  Number  Density 


VI .  Conclusion 


Sximmary  and  Conclusions 

The  results  of  this  study  show  that  in  an  inductively 
driven  plasma,  there  is  an  observable  phenomena  which  inter¬ 
acts  with  the  driving  Inductor  to  produce  changes  in  the 
electrical  circuit  impedance.  The  procedures  outlined  in 
Section  IV  can  be  used  to  identify  the  circuit  impedance 
which  exhibits  measurable  sensitivity  with  respect  to  changes 
in  gas  pressure,  input  power  or  driving  frequency.  Provided 
an  acceptable  circuit  model  of  the  plasma  is  introduced,  the 
measured  Impedance  can  be  used  as  a  diagnostic  to  determine 
the  plasma  parameters. 

In  pursuit  of  this  goal,  an  in  depth  theoretical  and 
experimental  analysis  was  made  of  a  specific  plasma  impedance 
measurement  technique  applied  to  an  inductively  excited  RF 
plasma.  Special  attention  was  given  to  calculation  of 
V' an,d  n.  Steady  state  plasma  theory  was  applied  to  explain 
predominate  electron  source  and  loss  mechanisms  and  basic 
plasma  phenomena. 

To  relate  the  plasma  to  the  electrical  driving  circuit, 
a  plasma-circuit  was  formulated.  Depending  heavily  on 
traditional  electric  circuit  analysis  but  recognizing  the 
nuances  of  RF  energy,  a  thorough  analysis  was  then  performed 
on  the  total  cii’cuit  as  well  as  individual  components.  Three 
results  evolved.  First,  the  method  of  data  reduction  out- 
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lined  by  Illc  (Illc«  D.B.,  1981:1544)  could  not  be  applied 
to  the  exiting  equipment.  Recognition  of  measurement 
sensitivity  and  error  propagation  indicated  that  inordinate 
amounts  of  RF  shielding  would  be  necessary  to  carry  out 
required  matching  network  calibration.  Second,  the  impedance 
measurement  procedure  described  in  Section  IV  was  conceived. 
This  procedure  offers  certain  advantages  over  that  presented 
by  Ilic.  Third,  equations  were  developed  to  calculate  plasma 
parameters  from  circuit  impedance  values.  For  comparison, 
the  suggested  development  adapted  in  Appendix  B  was  applied 
and,  although  some  comparable  data  was  obtained,  observable 
discrepancies  could  not  be  solely  attributed  to  either 
development  because  of  inherent  component  calibration  errors. 

More  than  Just  a  nvimerical  comparison  of  these  two 
methods  is  warranted  here.  A  subtle  difference  exists  which 
accents  each's  attributes  and  highlights  potential  applica¬ 
tions.  The  technique  described  in  Section  IV  utilizes  a 
direct  impedance  measurement  employing  the  RF  meter  described 
in  Section  III.  The  technique  presented  by  Ilic  and  adapted 
to  this  study  in  Appendix  B  requires  precalibration  of  matching 
network  components  in  order  that  control  knob  settings  may 
be  used  to  provide  data  input.  Both  techniques  presuppose  the 
Impedance  matched  condition  for  validity. 

It  is  possible  to  use  these  techniques  concurrently  (al¬ 
though  not  simultaneously)  so  as  to  make  two  measurements 
for  each  matched  Impedance  condition.  Illc's  method  is  per- 
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formed  on  an  operating  plasma  and  has  two  notable  advantages. 
It  provides  a  real-time  measurement  and  this  information  can 
be  used  as  an  error  signal  for  plasma  control  applications. 

The  method  conceived  in  this  study  requires  a  power  inter¬ 
ruption  but  offers  the  advantage  of  using  a  direct  circuit 
impedance  measurement  extracted  from  the  driving  inductor 
input  close  to  the  plasma,  thus  elminating  upstream  equipment 
calibrations.  Also,  being  a  post  plasma  measurement,  common 
stray  RF  propagation  effects  are  eliminated. 

Although  both  methods  were  applied  to  the  configuration 
under  study,  Ilic's  equipment  incorporated  a  continuous  tuning 
feature  which  would  destabilize  during  power  interruptions 
and  void  post  plasma  measurements. 

The  central  value  of  impedance  methods  is  that  they  are 
external  to  the  plasma  chamber  and  nondisturbing  to  either 
dynamic  particle  flows  or  temporal  and  spatial  fields.  This 
is  in  contrast  to  techniques  such  as  probes  or  microwave 
diagnostics  for  which  substantial  study  has  been  done  to 
account  for  these  effects.  Also  the  impedance  matched  cri¬ 
teria  imposed  assures  maximum  power  transfer  and  efficiency 
which  are  of  significant  benefit  to  industrial  users.  The 
contrasting  features  of  the  two  methods  offer  alternatives 
to  readers  concerned  with  specific  plasma  applications. 

The  plasma-circuit  model  Investigated  did  not  completely 
account  for  all  phenomena  occurring  and  does  not  permit  accu¬ 
rate  calculation  of  plasma  electron  number  density  or  electron 
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neutral  collision  frequency.  Associated  changes  in  stray 
capacitance  and  sheaths  impose  larger  effects  upon  the  cir¬ 
cuit  than  do  the  small  changes  in  plasma  parameters  of 
Interest . 

The  portion  of  the  model  presenting  the  most  difficulty 
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for  analysis  is  the  basic  capacitance,  C^.  The  magnitude 
of  this  term  is  small  (0.5  pf  in  this  experiment)  and 
difficult  to  obtain.  Furthermore,  changes  to  the  plasma 
result  in  only  a  percentage  change  in  (say  lOX).  The 
term  C^,  in  turn,  is  only  a  small  portion  (5)1)  of  the  total 
parasitic  capacitance  of  the  driving  inductor  and  it  is 
this  total  capacitance  that  interacts  with  the  electrical 
circuit.  This  interaction  itself  may  be  only  small  pertur¬ 
bation  at  frequencies  other  than  near  resonance.  The  result 
is  that  the  overa) 1  perturbation  of  the  circuit  impedance 
is  much  less  pronounced  than  the  real  changes  occurring 
within  the  plasma.  This  measurement  process  is  analogous  to 
determining  the  speed  of  a  car  by  measuring  the  Lorentz 
contraction  of  its  length. 

Recognizing  the  concomitant  demand  placed  on  measure¬ 
ment  and  calculation  accuracy,  a  detailed  error  analysis 
was  performed  on  the  primary  equations  used  to  calculate 
plasma  parameters.  The  equations  show  a  sensitivity  to 
errors  found  in  the  value  for  circuit  reactance.  Small 
measurement  errors  propagate  with  little  influence  to  the 
value  of  collision  frequency  and  a  lOX  confidence  level  may 
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be  assigned.  Calculated  number  density  on  the  other  hand 
is  strongly  influenced  by  errors  in  either  the  driving 
inductance  or  measured  reactance.  The  value  for  driving 
inductance  can  be  determined  with  high  accuracy  through 

successive  measurements,  however,  each  value  of  reactance 

« 

depends  upon  an  individual  circuit  impedance  value.  This 
sensitivity  allows  for  a  confidence  level  of  only  25%  for 
electron  number  density. 

In  addition  to  the  unresponsive  measurement  problem 
Just  mentioned,  there  exists  a  possible  model  defect  v^ich 
serves  to  further  distort  results.  The  total  model  capaci¬ 
tance  was  considered  to  consist  of  two  terms;  C  ,  the 

o 

capacitance  interior  to  the  tube,  and  Cg,  the  remaining 
parasitic  capacitance  of  the  system.  Both  terms  are 
determined  by  the  basic  geometry  of  the  system  and  are 
therefore  invarient  except  for  the  changes  in  already 
attributed  to  the  permittivity  function.  If  the  plasma 
forms  a  conducting  sheath  along  the  container  wall,  this 
may  introduce  an  additional  capacitive  term  between  the  in¬ 
ductor  and  the  plasma.  This  effect  was  noted  by  Robert  Olsen 
while  studying  the  effect  on  the  impedance  of  a  solonoid  due 
to  the  presence  of  a  center  conductor.  (Olsen,  R.A.,  1961:3) 
The  existence  of  such  a  term,  conductivity  dependant,  would 
vary  with  the  plasma  number  density.  Since  it  is  not  deter¬ 
mined  solely  by  geometry,  its  magnitude  cannot  be  measured 
by  the  same  technique  used  to  determine  the  terms  and  Cg. 
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Recommendations 


Tlte  procedures,  models,  and  results  developed  and  dis¬ 
cussed  in  this  study  coupled  with  the  preponderance  of  int¬ 
erest  shown  by  the  plasma  physics  community  warrants  further 
research  into  the  study, and  application  of  plasma  impedance 
diagnostics.  Research  should  include  both  capacitive  and 
inductive  reactors.  Measurement  sensitivity  and  the  dis¬ 
turbing  effect  of  stray  electrical  coupling  dictate  that 
inductively  oriented  research  be  slanted  toward  specific 
applications  where  parameters  such  as  frequency  or  power  can 
be  constrained  within  narrow  limits. 

In  addition,  the  simple  plasma-circuit  model  tested 
in  this  study  was  found  to  be  incomplete.  A  model  which 
would  better  account  for  observed  effects  is  shown  in 
Figure  22.  This  model  is  equivalent  to  the  model  tested 
over  a  limited  frequency  range  (Hayt,  W.H. ,  1962:458).  Also, 
the  resistor  in  series  with  the  inductor  would  account  for 
real  resistance  in  the  system  while  the  resistor  in  parallel 
with  the  capacitor  would  absorb  the  plasma  induced  resistive 
term. 

^1  I  C 

'  *'  —  o 
e _ 

.  T 


Figure  22.  Modified  Model 
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The  substitution  of  a  capacitive  type  plasma  reactor 
into  the  present  system  greatly  simplifies  the  model  needed 
and  offers  good  promise  of  producing  significant  results. 

It  is  likely  that  once  configured,  numerical  relations  could 
be  quickly  verified  using  a  Helium  plasma  before  continuing 
with  a  spectroscopic  study  of  a  molecular  gas.  The  design 
offered  in  Figure  23  is  currently  being  fabricated  by  the 
sponsor  of  this  thesis.  It  is  a  high  vacuum,  parallel  plate 
reactor  incorporating  a  micrometer  controlled,  5  cm  travel 
plasma  cavity.  Opposing  optical  windows  are  included  to 
facilitate  spectroscopic  analysis. 
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Appendix  A 


Data  Reduction 


The  equations  developed  here  represent  an  independant 
analysis  of  the  electrical  circuit  and  the  relationship 
between  measurements  and  plasma  parameters.  The  circuits 
shown  in  Figure  4  are  applicable  here.  We  recall: 


'P  ^  Vo^ 


(38) 


We  chose  to  write  the  permittivity  expression  in  the 
following  form: 


e^  *  1  - 
P 


v‘57 


(39) 


Combining  the  parallel  circuit  elements  and  making  the 
appropriate  substitutions  gives  the  expressions: 


vn  C 
e  o 


e  o 

+  V' 


(40) 


(41) 


The  term  ^  can  be  set  equal  to  some  real  and  imaginary  parts, 
e 

R*  and  X' ,  where  the  primes  indicate  that  the  values  relate 
e  e 

to  a  reciprocal  impedance.  Separating  the  parts,  solving 
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simultaneously  and  inserting  the  previously  defined  term  y 
produces  the  following  set  of  equations: 


V 


R'(u*  +  V*) 

■>.*  -Se; - 


(42) 


(43) 


values  for  R*  and  X'  are  readily  determined  from  a  measure- 
ment  of  and  the  use  of  a  calculator  program  labeled  "NUE" 
which  can  be  found  in  Appendix  C. 
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Appenctix  B 


Data  Redact log 


Tbe  original  concept  of  this  study  was  to  develop  a 
method  of  determining  plasma  parameters  from  the  matching 
network  dial  settings.  'This  was  accomplished  by  calibrating 
the  frequency  dependant  impedance  of  the  tuner  and  trans¬ 
mission  cable  and  borrowing  equations  developed  in  previous 
portions  of  this  study.  When  all  impedances  are  matched, 
the  signal  generator  drives  a  purely  resistive  load  of  * 
50  ohms.  Labeling  the  tuner  impedance  as  and  the  trans¬ 
mission  line  impedance  as  Z^t  we  can  write  the  equivalent 
impedance  as: 


Z^  »  Z^  -  Z.  -  Z„ 
e  o  t  m 


(44) 


Separating  into  real  and  Imaginary  parts  gives  the  following; 


50  -  (R^  +  R„) 
t  m. 


-  <*t  *  *m^ 


(45) 


(46) 


Once  the  desired  frequency  is  selected,  a  measurement  of  the 
transmission  line  gives  a  value  for  the  associated  impedance. 
Extensive  use  of  this  method  would  suggest  calibrating  the 
transmission  line  over  a  continuous  range  of  frequency  there¬ 
by  avoiding  replicating  '.his  measurement.  The  matching  net- 
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work  dial  settings  are  read  at  the  impedance  matched  point. 
These  dial  settings  correspond  to  values  of  the  tuner  com¬ 
ponents,  C^,  C2*  ^  obtained  from  calibration  curves 

shown  in  Figures  24,  25,  and  26.  These  values  can  be  re¬ 
duced  by  circuit  analysis  (a  short  computer  program  was 
written  to  perform  repeated  calculations)  to  yield  a  value 
for  the  combined  tuner  impedance.  Equations  (45)  and  (46) 
then  give  values  which  reduce  to  a  single  impedance  given 
by  the  equation: 

-uLX  -  IcdLR  . 

7  a  _ ^  ®  (47^ 

P  -Ro  + 

The  real  and  imaginary  parts  of  are  then  used  in  the 
calculator  program  labeled  "ZIP'*  (See  Appendix  C)  to  calcu¬ 
late  values  for  collision  frequency  and  number  density. 


400  h  A  30  MHZ 


Figure  25.  Capacitor  Calibration 


Appendix  C 
Computer  Programs 


The  following  two  program  are  designed  for  use  on  the 
Hewlett-Packard  Model  HP-41C  programmable  calculator. 

Program  "NUE"  or  "ZIP*'  ‘ 

Initialize  program 

Store  w  in  register  06. 

Store  the  phase  of  In  register  Y. 

Store  the  magnitude  of  in  register  X. 


Program : 

01. 

LBL  NUE 

02. 

X— y 

• 

CO 

o 

CHS 

04. 

X— Y 

05. 

P-R 

06. 

XROM  CINV 

07. 

STO  11 

08. 

RDN 

09. 

STO  12 

10. 

RCL  06 

11. 

X" 

12. 

STO  13 

13. 

2.70  E-6 

14. 

* 

15. 

1/x 

16. 

CHS 

17. 

8.92  E-12 

18.  . 

+ 

72 


19. 

RCL  06 

20. 

* 

21. 

CHS 

22. 

RCL  12 

23. 

+ 

24. 

1/x 

25. 

RCL  11 

26. 

27. 

RCL  06 

28. 

* 

29. 

STO  02 

30. 

NUE 

31. 

ARCL  X 

32. 

AVIEW 

33. 

RCL  02 

34. 

35. 

RCL  13 

36. 

+ 

37. 

RCL  11 

38. 

* 

39. 

STO  14 

40. 

3.182  E9 

41. 

RCL  02 

42. 

* 

43. 

0.52  E-12 

44. 

45. 

RCL  14 

46. 

X— Y 

47. 

/ 

48. 

STO  03 

49. 

Ne  = 

50. 

ARCL  X 

51. 

AVIEW 

52. 

END 

I 
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Proeram: 

01. 

LBL 

ZIP 

02. 

X—Y 

03. 

CHS 

04. 

X—Y 

r 

05. 

P-R 

06. 

STO 

11 

07. 

RDN 

08. 

STO 

12 

09. 

Store  Rx 

10. 

PROMPT 

11. 

RCL 

06 

12. 

2.70  E-6 

13. 

* 

14. 

STO 

13 

15. 

RCL 

12 

16. 

- 

17. 

STO 

14 

18. 

RCL 

12 

19. 

RCL 

13 

20. 

21. 

CHS 

22. 

STO 

15 

23. 

RCL 

11 

24. 

RCL 

13 

25. 

* 

26. 

RCL 

15 

27. 

RCL 

14 

28. 

RCL 

11 

29. 

CHS 

30. 

XROH  C/ 

31. 

STO 

11 

32. 

RDN 

33. 

STO 

12 
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34. 

17.2 

35. 

STO 

13 

36. 

0.52  E-12 

37. 

STO 

14 

38. 

RCL 

06 

39. 

* 

40. 

STO 

17 

41. 

1/x 

42. 

RCL 

12 

43. 

* 

44. 

STO 

15 

45. 

RCL 

12 

46. 

47. 

RCL 

13 

48. 

* 

49. 

STO 

16 

50. 

RCL 

11 

51. 

X* 

52. 

RCL 

13 

53. 

* 

54. 

RCL 

16 

55. 

+ 

56. 

RCL 

15 

57. 

+ 

58. 

STO 

15 

59. 

RCL 

14 

60. 

* 

61. 

1/x 

62. 

RCL 

11 

63. 

« 

64. 

STO 

07 

65. 

NUE 

« 

66. 

ARCL  X 

67. 

AVIEW 

68. 

RCL 

12 
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RCL  13 


RCL  17 


STO  18 
RCL  17 
RCL  15 

* 

1/x 

RCL  18 
* 

RCL  13 
+ 

RCL  06 


3.182  E9 


Ne  = 
ARCL  X 
AVIEW 


m. 
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